Partial wave analysis of the reaction 7p poo and the search for nucleon resonances 
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An event-based partial wave analysis (PWA) of the reaction 7p — > puj has been performed on a 
high-statistics dataset obtained using the CLAS at Jefferson Lab for center-of-mass energies from 
threshold up to 2.4 GeV. This analysis benefits from access to the world's first high-precision spin 
density matrix element measurements, available to the event-based PWA through the decay distri- 
bution ol uj ^ ■K^-K~-K^ . The data confirm the dominance of the t-channel tt" exchange amplitude 
in the forward direction. The dominant resonance contributions are consistent with the previously 
identified states _Fi5(i680) and Z)i3(1700) near threshold, as well as the Gi7(2190) at higher ener- 
gies. Suggestive evidence for the presence of a = 5/2^ state around 2 GeV, a "missing" state, 
has also been found. Evidence for other states is inconclusive. 

PACS numbers: 11.80.Cr,11.80.Et,13.30.Eg,14.20.Gk,25.20.Lj,23.75.Dw 



I. INTRODUCTION 



Studying near-threshold lo photoproduction presents 
an interesting opportunity to search for new baryon res- 
onances. Measurements made by previous experiments 
have produced relatively high-precision cross sections at 
most production angles; however, precise spin density 
matrix elements have only been measured at very for- 
ward angles and at higher energies [TJ |2J El HI [5] . In the 
near-threshold region, the only previously published spin 
density matrix results, which come from the SAPHIR 
collaboration, constitute a total of 8 data points in the 
energy range from oj photoproduction threshold up to a 
center-of-mass energy, W, of about 2.4 GeV [5j. 

A number of theoretical studies have been undertaken 
with the goal of extracting resonance contributions to 
LO photoproduction from these data. All of the authors 
agree on the importance of contributions from tt*' ex- 
change in the f-channel; however, discrepancies exist on 
the importance of various resonance contributions. In 
the calculations of Oh et al. [B], the dominant reso- 
nance contributions are found to be from a "missing" 
Pi3(1910) state (i.e. a state predicted by the constituent 
quark model but not observed experimentally) and from 
a Z3i3(1960) state. In contrast to this, Titov and Lee 
find the most significant resonance contributions to be 
from the 1513(1520) and Fi5(1680) states. The quark 
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model calculations made by Zhao [5] find that the two 
most important resonance contributions to to photopro- 
duction come from the Pi3(1720) and Fi5(1680) states. 
The Pii(1710) and Pi3(1900) states were found to be the 
dominant resonance contributions in the coupled-channel 
analysis of Penner and Mosel [S] . 

All of the models mentioned above were fit solely to 
differential cross sections. A more recent analysis [TD] 
also included the spin density matrix elements published 
by the SAPHIR collaboration This work found 

the largest resonant contributions to oj photoproduction 
to be from the sub-threshold 1515(1675) and Fi5(1680) 
states. The authors noted the importance of the strong 
additional constraints placed on their model by the polar- 
ization information and concluded that: there is urgent 
need for precise measurements of the spin density ma- 
trix in more narrow energy bins to pin down the reaction 
picture. 

Recently published results from the CEBAF Large Ac- 
ceptance Spectrometer (CLAS) have provided such mea- 
surements [IT]. In the center-of-mass (cm.) energy range 
from threshold up to 2.84 GeV, differential cross section 
results were reported at 1960 points in W and cos0"„ . 
The experiment did not use a polarized beam or a polar- 
ized target; thus, only the Pqq, Pi-i and i?e(pio) elements 
of the spin density matrix could be determined (the defi- 
nitions of which can be found in |12j). These results were 
reported at 2015 points in W and cosO'^j^ . The increase 
in precision for p\jm'-< energy range overlapping 

the SAPHIR results, is approximately a factor of 148. 

In this paper, we present an event-based mass- 
independent partial wave analysis (PWA) of these data, 
i.e. the data are only divided into narrow cm. energy 
bins. In each of these narrow bins, the spin-independent 
part of any resonance propagator - a complex function 
of - is approximated as a constant complex number. 
This allows us to reduce model dependence in our treat- 



ment of resonances (see Section IIIC for further discus- 
sion on this topic). 

The data used in our analysis were obtained using the 
CLAS housed in Hall B at the Thomas Jefferson Na- 
tional Accelerator Facility. Real photons were produced 
via bremsstrahlung from a 4.02 GeV electron beam. The 



3 



momenta of the recoihng electrons were then analyzed in 
order to obtain the energy of the photons with an uncer- 
tainty of 0.1% [T^. The physics target was filled with liq- 



uid hydrogen. The ui 



decay was used to select 



the reaction of interest. The momenta of the charged par- 
ticles (p, TT"*" , 7r~ ) were determined using the CLAS detec- 
tor with an uncertainty of approximately 0.5%. The neu- 
tral tt" was reconstructed using kinematic fitting. More 
details concerning the analysis techniques can be found 
in [TTJ [Tl] . A detailed description of the CLAS can be 
found in [15] . 

In total, the dataset consists of over 10 million signal 
events divided into 112 10-MeV wide cm. energy bins. 
Our primary interest is in extracting possible nucleon res- 
onance contributions; thus, we have restricted our PWA 
to include only bins with W < 2.4 GeV. In total, 67 cm. 
energy bins were used in the PWA (the W = 1.955 GeV 
bin was excluded due to issues with the normalization cal- 
culation [TT ) . This work represents the first event-based 
PWA results on baryons from photoproduction data. 



II. PWA FORMULAS 

As stated above, to limit theoretical model dependence 
we divided our data into 10-MeV wide W bins. Thus, 
all formulas written below are intended to describe data 
from a narrow cm. energy range. In all of the work that 
follows, pi,pf,q and k will be used for the initial proton, 
final proton, iv and photon 4-momenta. The z-axis in the 
overall cm. frame, defined by fc, is used as the angular 
momentum quantization axis. The Mandelstam variables 
are defined as 



s = {p, + kf = {pf + qf 

t = {q-k)' = {p.,-pff 

2 



U 



(la) 
(lb) 
(Ic) 



The mass of the proton and lo are denoted as Wp and Wu} , 
respectively. 

The Lorentz invariant transition amplitude, Ai, of the 
process 7p puj — )■ pTr+vr^Tr*^, can be written as 



|7W(a, x)f 



E 



rrii ,m~. .rrif 



^ ^ '^mi, 771^ ,7nf (^7 *^) 



(2) 



where mi,m~f,mf are the initial proton, incident pho- 
ton and final proton spin projections, A'l^y^. are the 
partial wave amplitudes (the form of which is discussed 
in Section III I , x denotes the complete set of kinematic 



variables describing the reaction and a. are the unknown 
parameters to be determined by the fit. We denote the 
detector acceptance by rj{x) and the phase-space volume 
as d^{x) = (j){x)dx. A more detailed description of the 
work presented in this section is given in 14J. 



A. Likelihood 

All of our fits are event-based; thus, the data were only 
binned in W. To obtain estimators for the unknown pa- 
rameters, a, we employ the extended unbinned maximum 
likelihood method. The work detailed in this section is 
based on that of Chung [TB]; however, the normalizations 
we have developed differ from his work. The likelihood 
function is defined as 



£ = 



]^7'(a,fj), 



(3) 



where the term in parentheses is the Poisson probability 
of obtaining n events given the expected number h(d) 
(the calculation of which is discussed below), Xi rep- 
resents the complete set of kinematic variables of the 
j*'' event and V{d,x) is the probabihty density function 
given by 



V{d,Xi) 



\M{d,Xi)\'^if]{xi)(l){xi) 
J \M{d, x)\'^ri{x)4){x)dx' 



(4) 



^From left to right, Eq. Q accounts for the relative 
strength of the transition amplitude, the detection prob- 
ability and the available phase space for the i*^ event. 
Calculation of the denominator, which normalizes the 
probability density function, is discussed below. The es- 
timators d are then found by maximizing C. 



B. Normalization 

The expected number of signal events for a given set 
of parameters is given by 

n{d) = / \M{d,x)\-'7j{x)d^x), (5) 



8(s - wl) 

which includes the average over initial spin states. 



T{s) 



^{s)ptarg£targNAb 



^targ 



(6) 



is the "target factor" obtained from the target density, 
Ptarg, length, itarg, and atomic number, Atarg] along 
with Avogadro's number, iV^; the branching fraction of 
Lu —)■ tt+tt^tt'', b; and the integrated photon flux in each 
W bin, JP(s). 

The integral in Eq. ([5| must be done numerically due 
to the lack of an analytic expression for the detector ac- 
ceptance. Monte Carlo events were generated in each W 
bin according to jp — > puj (the lo mass was generated 
according to a Breit-Wigner distribution), uj —> 7r"'"7r~7r'^ 
phase space and then run through a GEANT-based de- 
tector simulation package (discussed in detail in [TTJ[T3]). 
This procedure simulates the acceptance of the detector 
by rejecting events that would not have survived the data 
analysis, i.e. for each generated event, the acceptance 
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factor ri(xi) = or 1. The integral can then be approxi- 
mated by 



\M{d, x)\^ri{x)d^{x) 



l^j:\Mia,x.)\^ (7) 



where ngen{nacc) is the number of generated (accepted) 
Monte Carlo events and 



d^{x) 



[{s - {Wp + w^)^){s - (Wp - w^)^)] 
4(27r)5s 



(8) 

is the volume of the 2-body puj phase space (the Stt phase- 
space volume is factored into the normalization of the u 
decay amplitude). 

Using Eqs. ([t]) and Eq. ^ can be rewritten as 



n{d) 



[{s - (Wp + w^)^){s - (Wp - w^)^)] 



1/2 



647rs(s — w^) 

r(.)'^ 



''gen 



J2\Mid,x-.)\\ (9) 



This normalization allows us to use physical coupling 
constants in our event-based fits, i.e. it allows us to put 
our parameters on an absolute scale. Thus, our normal- 
ization scheme permits direct theoretical input. 



C. Log likelihood 

Due to the monotonically increasing nature of the nat- 
ural logarithm, the likelihood, C, defined in Eq. ^ can 
be maximized by minimizing 

n 

— ln£ = — nln h{d)+lnnl + fi{d) In 7^(0?, x^), (10) 

i 

which, using Eqs. (|4]) and ([5|, can be rewritten as 

8(s-u)2) 

- In/: = Inn! + n{d) + nln 

T(s)(27r)4 

n 

-^ln|A^(a,x,)pr7(x,)0(x,). (11) 



D. Handling background 

To accurately extract partial wave contributions to u 
photoproduction, background events, i.e. non-w events, 
must be separated from the signal in a way that pre- 
serves all kinematic correlations. The method we applied 
to our data, described in detail in [TTl [TH |T7] , assigned 
each event a signal quality factor, or Q-factor. This back- 
ground was assumed to be non-interfering. Following our 
previous work |17j . we can rewrite Eq. (12) using these 
Q-factors as 

n 

- ln£ = - ^ gjn \M{d,Xi)\'^ + n{d) + const, (13) 



where Qi is the Q-factor for the i event. Thus, the 
Q-factors are used to weight each event's contribution to 
the likelihood. We also note that in the literature, the t- 
and M-channel contributions are often referred to as back- 
ground. This theoretical background is not to be confused 
with the experimental background discussed here. In this 
analysis, the theoretical backgrounds were allowed to in- 
terfere with the s-channel amplitudes in our PWA. 



III. PWA AMPLITUDES 

The choice of which amplitudes to include to describe 
the data is partially motivated by experimental measure- 
ments. The Lu photoproduction cross section is known to 
have a strong forward peak, even at near-threshold ener- 
gies [IT]. At higher energies, the cross section develops 
a rather pronounced backwards peak as well |llj . These 
features are typically associated with meson and nucleon 
exchanges in the t- and u-channel, respectively. The re- 
cent CLAS data also possess a number of features in the 
cross sections and spin density matrix elements sugges- 
tive of resonance contributions |llj . Thus, it would seem 
that S-, t- and u-channel amplitudes may be required 
to fully describe the data (see Fig. [T]). The formalism 
used to construct these amplitudes is described in de- 
tail elsewhere |14j . below we simply give an overview of 
the different types of amplitudes used in our analysis. 
All of these amplitudes were computed using the qf t++ 
package [18j , which performs numerical computations of 
quantum field theory expressions. 



Neglecting terms that do not depend on the parameters, 
we can then rewrite Eq. (Ill as follows: 



-InC^ - ^ln\M{d,Xi)\'^ + n{d) + const. (12) 

i 

We note here that for any set of estimators that minimize 



Eq. (12), the expected number of events is n(d) — n. 



The Lo 



amplitude can be written in terms 
of the isovectors, and the 4-momenta, p^r, of the pions, 
along with the uj 4- momentum (q), polarization (e), and 
spin projection, m^^ as 



xep^a.f}p'^+p'^^p'l„e''{q,m^), (14) 
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FIG. 1; Feynman diagrams for the amplitudes used in our 
analysis. The images were produced using the JaxoDraw 
package [19]. 



which is fully symmetric under interchange of the three 
pions. For this reaction, where all final states contain 
uj — > 7r"'"7r~7r'', the isovector triple product simply con- 
tributes a factor to the global phase of all amplitudes. In 
the Lo rest frame, Eq. (14) simplifies to 



-^"-^tt+tt-ttO « iP^+ X pV-) • e(m^), (15) 
which is the standard non-relativistic result PU] , 

B. t- and tt-channel 

Previous studies of forward lu photoproduction data 
have shown that the reaction is dominated at low energies 
by pion exchange and at higher energies by diffractive 
processes, i.e. Pomeron exchange (see, e.g., [T]). We 
have chosen to use the non-resonant terms included in 
the model of Oh, Titov and Lee (OTL) [H] in our partial 
wave analysis. The OTL Pomeron exchange amplitude 
follows the work of Donnachie and Landshoff [21 with 
the unknown parameters fixed by fitting to high energy 
vector meson cross section data. 

The OTL model also includes pseudoscalar meson 
exchange amplitudes obtained from the following La- 
grangians: 



(16a) 
(16b) 



where (/) = (vTjTy), and ^ denote the pseudoscalar, 
photon and proton fields, respectively. The vertices in 
these amplitudes were dressed using form factors of the 



F(t,A) = 



A2 



W0 



(17) 



where A is the cutoff parameter for the interaction and 
is the mass of the exchanged particle. The larger 
mass and weaker coupling constants of the r] suppress its 
contribution relative to that of the pion. 

Nucleon pole terms were obtained from the La- 
grangians 



2m 



(18a) 



LUpp 



= -gupp^h"" - ^a^'^d,)^^^. (18b) 



The form factor 



Air 



A^ 



(19) 



where x = {s,u), was included to dress the corresponding 
vertices. The amplitudes were also modified to preserve 
gauge invariance. The details concerning these modifi- 
cations, along with the values of the parameters used in 
the model can be found in [B] . 

In the near-threshold region, the high precision spin 
density matrix results published by CLAS confirm the 
dominance of ^-channel 7r° exchange in the forward di- 
rection; however, at higher energies the existing theoret- 
ical models fail to reproduce the CLAS data [22j. The 
Pomeron amplitudes are able to describe the energy de- 
pendence of the forward cross section, but fail to ade- 
quately describe the spin density matrix elements. The 
unknown parameters present in the nucleon exchange 
amplitudes can be modified to describe the backward- 
angle data at higher energies if some assumptions, the 
reliability of which are unknown, are made |22j . 

Our analysis is restricted to the energy range from 
threshold up to 2.4 GeV. For cm. energies below 2 GeV, 
the tt" exchange amplitude dominates the t-channel con- 
tributions. In the higher energy range used in our 
analysis, the tt" contribution is still substantially larger 
than that of the Pomeron. Thus, the deficiencies in 
the Pomeron amplitudes (discussed above) should not 
greatly affect our PWA results. For this reason, we have 
chosen to use the OTL i-channel terms with the param- 
eter values obtained in that analysis [B]. Due to the un- 
reliability of the assumptions under which the nucleon 
exchange parameters were determined, we have decided 
not to include these amplitudes in our analysis, i.e. we do 
not include any u-channel terms in our fits. The effect of 
this choice on our conclusions was found to be negligible 
(see Section [v| . 

C. Resonant waves 

The formalism used to construct our resonant ampli- 
tudes is described fully in [14 . It involves the use of 
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relativistic tensor operators and is similar to the frame- 
work employed by Anisovich et al. |23j . 

As discussed above, we do not impose resonant-like 
shapes on our s-channel waves. Instead, we divide our 
data into narrow cm. energy bins. In each of these bins, 
the mass-dependence of an s-channel wave with spin J 
and parity P, which we will denote TZjp(s), is approxi- 
mated by a constant complex number: 



^r^pe^*.-e(<5-|yi-I¥fc|), (20) 

b 

where the sum is over the 10-MeV wide cm. energy 
bins, r''jp,(j)^jp are the strength and phase of the mass 
dependence in each bin, respectively, and <5 = 5 MeV is 
the maximum distance in any bin from the centroid Wb- 
The resonant waves then enter into our fits according to 

y^mi,m^,m; — 9 LSi9 LS f ''^ J'' \^ I -^"mi,™,^ ,m f , {■^^) 

where LSn^fj are the angular momentum quantum num- 
bers of the initial (final) state, are the unknown 
coupling constants to these states and A is the covariant 
ampUtude obtained using the formalism described in [14]. 

The values extracted for each s-channel wave's r jp and 
parameters can be used to search for evidence of nu- 
cleon resonance contributions in that wave. For every fit 
iteration run in each W bin (multiple iterations are run 
to alleviate problems caused by local minima), these pa- 
rameters are started at random values that include the 
entire physically allowed range of the parameter. For 
example, (f)^jp is started randomly in the range [0,27r). 
Estimators for the parameters are then found by maxi- 
mizing the likelihoods independently in each W bin. In 
this way, the mass dependence of the waves is extracted 
unbiasedly. If resonant-like features are found in the cross 
sections and phase motion of the s-channel waves, then 
this is very strong evidence that resonances do contribute 
to the scattering amplitude. 

If the strength observed in an s-channel wave is due 
to a single resonant state, then TZjp{s) should be (at 
least qualitatively) described by a const ant- width Breit- 
Wigner line shape of the form: 

= ?r-T' (22) 

where w and F denote the mass and width of the state, re- 
spectively. If, however, the strength in the wave is caused 
by multiple resonant states or from a non-resonant pro- 
cess, then the use of a Breit-Wigner line shape is not 
valid. The line shape given in Eq. ( 22 ) neglects the kine- 



matics and dynamics of the mass dependence of the res- 
onance. This can be an issue, e.g., near threshold. This 
will be addressed below. 



IV. RESULTS 

Before examining our results, it is important to re- 
iterate the goals of our analysis. We are attempting to 



extract strong resonance contributions to u photoproduc- 
tion in a least model-dependent way. We do not enforce 
resonance shapes on the mass dependence of our partial 
waves. Instead, we first extract the strengths and phase 
motion of our partial waves independently in each W bin. 
This stage will be referred to as a partial wave extraction 
(PWE). The second stage involves comparing the results 
of the PWE to what is expected from resonances in a 
mass-dependent fit (MDF). The PWE's are performed 
using unbinned extended maximum likelihood fits to the 
data in each W bin. The MDF's are simply fits to 
the phase differences obtained between partial waves in 
a given W region. 

We are not looking to build a complete model of oj pho- 
toproduction, i.e. we do not claim that the fits discussed 
below contain all of the amplitudes that contribute to 
this reaction. Because of this, we do not expect the phys- 
ical observables extracted from our fits to provide perfect 
descriptions of our measurements; however, the descrip- 
tions in many kinematic regions are very good. Finally, 
we are not attempting to extract all resonance contribu- 
tions to Lo photoproduction, only the most significant. 



A. Choice of wave sets 

It is important to have a systematic method for se- 
lecting wave sets. As discussed in Section IIIB all of 



the wave sets used in our analysis contain the OTL t- 
channel terms (which contain no free parameters) and 
no u-channel terms. Systematic studies show that the ef- 
fects on extracted resonance parameters due to the choice 
of the non-resonant model are small (see Section |v]) . 

The mass-independent nature of our procedure, i.e. 
the bin-to-bin freedom of the resonance parameters, 
makes the use of smaller wave sets advantageous. 
For this reason, we began our wave selection pro- 
cess by scanning the entire energy range of interest, 
1.72 GeV< W <2.4 GeV, using the OTL t-channel terms 
along with waves from a single spin-parity, . The 
goal of this scan was to identify (possible) energy ranges 
where waves of a given perform significantly better 
than waves of any other spin-parity. This information 
alone does not constitute evidence of resonance produc- 
tion; however, it can serve as a guide as to which waves 
are more likely to contribute strongly to uj photoproduc- 
tion. 

Given a pair of fits run with different wave sets, the 
difference in the log likelihoods obtained from the fits, 
A In £ = In £(j ~ In can be used to quantitatively de- 
termine which fit best describes the data. If Aln£ > 0, 
then wave set a provides a better description of the data 
than wave set b, while A In £ < implies the converse is 
true. 

Figure |2] shows two examples comparing the likelihood 
differences between fits with different s-channel waves. 
The A In £ quantities are shown for two separate fits, one 
with = 5/2+ and one with = 3/2^, each of which 
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is compared to a fit with — 1/2+ . ^From threshold up 
to - 1.85 GeV, the fit with = 3/2" is clearly the 
best, while in the energy range 1.85 GeV< W <2 GeV 
the preferred wave is — 5/2+. It is also clear in Fig. 2] 
that both the = 3/2^ and = 5/2+ waves provide 
better descriptions of our data than the — 1 /2+ wave 
in this energy range. 

Similar fits were run using any single s-channel wave 
with J < 9/2 of both parities. In the region from 
threshold up to - 1.85 GeV, the = 3/2" wave 
was found to provide a better description of our data 
than any other wave. Similarly, in the energy range 
1.85 GeV< W <2 GeV, the = 5/2+ wave was found 
to provide the best description. 

Scans were also performed using two s-channel waves 
and the OTL t-channel terms. The quantity Aln£ can 
also be used in these fits to determine which wave sets 
best describe our data. In the W <2 GeV region, the 
best fit was obtained using the s-channel waves with 
JP = 3/2", 5/2+. This wave set had the best hkeli- 
hood in every bin in this energy range (typically by a 
large amount). Given the results of the single wave s- 
channel scans discussed above, this is not a surprising 
result. Above 2 GeV, the preferred wave set consisted 
of the — 5/2+,7/2~ waves, along with the t-channel 
terms. As in the lower energy range, this wave set had 
the best likelihood in every energy bin for W > 2 GeV. 

The results presented for the waves below were not af- 
fected by our choice of wave set; however, as the number 
of waves was increased so did the noise. For this rea- 
son, we have chosen to present the results from fits with 
at most three s-channel waves. See Section |V] for more 
discussion on fits with a larger number of waves. 
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FIG. 2: Aln£ vs W^(MeV): Example likelihood differences 
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B. Fit I: the near-threshold region 

Our preliminary s-channel scans showed that 
the best fit using two s-channel waves, along 
with the OTL i-channel terms, in the energy 
range 1.72 GeV <W< 2 GeV is obtained using 
= 3/2- , 5/2+ . To extract any possible resonance 
contributions in these waves, event-based PWE fits 
were run using the locked OTL i-channel terms along 
with jP = 3/2~,5/2+ s-channel waves parametrized 
according to Eq. (21 1. In each energy bin, multiple fit 



iterations were run using random starting values for the 
parameters; the results presented below are always from 
the fit with the best likelihood. 



1. Cross sections and phase motion 

The strength and phase of each s-channel wave were 
completely free to vary in each energy bin, i.e. they were 
fit independently. The cross sections extracted for the 
s-channel waves are consistent with either near or sub- 
threshold resonance states (see Fig. [3|. The PDG lists 
two states in these waves consistent with this hypothe- 
sis: (1) the 4-star i^i5(1680), which has a well known very 
large coupling to 7p; (2) the 3-star I?i3(1700), which is 
currently rated as having only a 2-star coupling to 7p. 
We note here that the masses of the states cannot be 
(precisely) estimated by simply examining the cross sec- 
tions due to threshold suppression effects. 

Figure |4] shows the phase motion between the two s- 
channel waves obtained from the PWE fits. The phase 
differences were then fit in a MDF using the constant 
width Breit-Wigner line shapes of the form given in 
Eq. p2| ). We chose not to use mass-dependent widths 
in the Breit-Wigner line shapes (despite the proximity 
to puj threshold) since the J^i5(1680) (and perhaps the 
Di3(1700) as well) is below threshold, which introduces 
model dependence in a single-channel analysis. For this 
reason, extracting precise values for the resonance param- 
eters may not be possible; however, the use of Eq. (22 1 is 



sufficient to provide evidence for the presence of known 
PDG states in our data. 

In principle, our MDF fits could have been made to 
the partial-wave intensities as well as the phase differ- 
ence. In order to do this, form factors would need to 
be introduced at both the production and decay vertices. 
While these form factors do not strongly infiuence the 
phase difference, they are very important for describing 
the shapes of the cross sections. The intensities of the 
partial waves (in this fit, and in the following sections) 
are qualitatively consistent with the expected resonance 
shapes. Obtaining good quantitative agreement requires 
the extra degrees of freedom introduced by the form fac- 
tors; however, including these factors also introduces ad- 
ditional model dependence. Thus, we have decided to 
only fit the phase difference — which, as noted, is nearly 
independent of the form factors. 
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Our results are, qualitatively, in good agreement with 
those expected from the PDG states mentioned above. 
The dashed-line in Fig. |4] was fit requiring all param- 
eters to be within the limits quoted by the PDG for 
the Fi5(1680) and Di3(1700). There is a minor dis- 
crepancy in the near-threshold bins. The parameters 
of the Z3i3(1700) are not as well known as those of the 
Fi5(1680); thus, we also performed a MDF allowing the 
3/2~ parameters to vary freely. This fit resulted in a mass 
of 1754 MeV and a width of 39 MeV for the 1)13(1700), 
which are very close to the PDG limits. Uncertainties 
of 21 MeV for the mass and 12 MeV for the width were 
estimated by examining the variation in the 3/2~ pa- 
rameters while using various parameter values (all within 
the PDG limits) for the 5/2"'", along with fitting different 
sub-ranges in W of the phase motion. 
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FIG. 4: Results from Fit I: A<j!> — 03/2- — (^5/2+ (radians) vs 
W{MeV): The dashed line was fit using constant width Breit- 
Wigner distributions requiring the parameters to be within 
the limits quoted by the PDG for the 4-star _Fi5(1680) and 
3-star Di3(1700). The solid line was fit allowing the 3/2~ 
parameters to vary freely, the results are listed in the text. 
The error bars are purely statistical. 
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FIG. 3: (Color Online) Results from Fit I: a{fj.b) vs VF(MeV): 
Total cross sections from all of the waves included in the fit 
(filled squares), only t-channel waves (open squares), only 
— 5/2^ waves (circles) and only — 3/2" waves (trian- 
gles). The cross sections extracted for both s-channel waves 
are consistent with near/sub-threshold resonances. The error 
bars are purely statistical. 



A single channel analysis is not the best environment 
for extracting precise resonance parameters; however, the 
qualitative agreement of the phase motion obtained from 
the PWE's to that of the two PDG states is very sug- 
gestive of their presence in our data. Recall that the 
phase parameters were each started pseudo-randomly in 
the range [0, 2tt) in each W bin. Yet the results are in 
good qualitative agreement with the phase motion ex- 
pected using simple constant width Breit-Wigner distri- 
butions for the PDG i^i5(1680) and 1)13(1700) states. 



2. Production helicity amplitudes 

We can also compare the production helicity couplings 
extracted from our fits to those quoted by the PDG. The 



ratio of the helicity amplitudes is obtained using the pro- 
duction couplings, gj^g. in Eq. (21), extracted by the fits 
and the s-channel production amplitudes. Due to the 
nature of the covariant formalism, these amplitudes are 
energy-dependent; i.e. the ratio of the helicity ampli- 
tudes is a function of W . The PDG quotes these values 
at the resonance masses. For the Di3(1700), the PDG 
reports the ratio of the helicity amphtudes as [53] 



A 



3/2 



Al/2 



= 0.11 ± 1.34. 



(23) 



Our fits extract this value to be in the range [—0.06, 0.13], 
depending on the mass of the = 3/2~ state. This is 
in a very good agreement with the PDG value. 

The ratio of the helicity amplitudes for the Fi5(1680) 
extracted by our fits is consistent with the PDG [23] 
value. However, projecting this ratio from the pto thresh- 
old down to the required mass makes obtaining a precise 
quantitative value difficult. 



3. Comparison to observables 

Fit I consists of three production mechanisms: 

(1) OTL t-channel terms (with no free parameters), 
which are dominated by ir'^ exchange in this energy range; 

(2) — 3/2~ s-channel waves, whose extracted pa- 
rameters are consistent with the PDG Di3(1700) state; 

(3) = 5/2"*" s-channel waves, whose extracted param- 
eters are consistent with the PDG i^i5(1680) state (at 
least, near threshold). This is almost certainly not all 
of the physics contributing to lo photoproduction in this 
energy range. Thus, we would not expect Fit I to provide 
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a perfect description of our data. Before we examine the 
quahty of Fit I, we note that the OTL i-channel terms 
do provide a good description of our forward data in the 
energy regime of Fit I. 

Figures |5] and [6] show comparisons of the differential 
cross sections and spin density matrix elements extracted 
from the PWA fits compared to our measurements [TT]. 
Recall that we do not fit to the experimental observables 
directly. We perform event-based fits to the data used to 
obtain the measured results. The forward cross section 
and polarization observables are very well described in 
this energy range, confirming that 7r° exchange in the 
t-channel does dominate the amplitude in this region. 
There is a discrepancy in the description of the cross 
section at backwards angles that increases with energy. 
This could be due to the lack of it-channel terms. It 
could also be a signature of unaccounted-for s-channel 
amplitudes. 

Though we did not start off by including known PDG 
resonance states, the fit has extracted evidence for them 
from our data. We also note here that the large-angle 
cross section at W — 1.8 GeV is virtually flat. With- 
out polarization information, the production mechanism 
could have easily been mistaken for a J = 1/2 wave. This 
demonstrates the importance of the spin density matrix 
elements. 

The quality of the description of the observables de- 
creases slightly with increasing energy. This signifies that 
there is another production mechanism that is not ac- 
counted for in the fit . This is expected due to the limited 
number of waves included in the PWA. Adding additional 
waves improves the description of the data but has vir- 
tually no effect on the strengths and phase motion of 
the two s-channel waves presented in this section. Thus, 
the conclusions drawn about resonance contributions are 
robust, and do not change when additions are made to 
the wave set. See Section [IV D 1| for results obtained by 
adding an additional s-channel wave to this fit. 



C. Fit II: the higher mass region 

Our preliminary s-channel scans showed that the best 
fit using two s-channel waves and the OTL f-channel 
terms in the energy range 2 GeV< W < 2.4 GeV is ob- 
tained using waves with = 5/2+, 7/2^. To extract 
any possible resonance contributions in these waves, the 
same procedure used for Fit I was employed. 



1. Cross sections and phase motion 

As in Fit I, the strength and phase of each s-channel 
wave were completely free to vary in each energy bin 
in the PWE's. The cross sections extracted for the s- 
channel waves, shown in Fig. [7] are consistent with the 
tail of a = 5/2+ state with a mass below 2 GeV (as 
seen in Fit I) and a = 7/2~ state with a mass near 



2.2 GeV. The PDG lists a state consistent with this hy- 
pothesis: the 4-star Gi7(2190), which currently is only 
listed as having 1-star coupling to jp. 

Figure [8] shows the phase motion between the two s- 
channel waves extracted from the PWE's. Our results do 
not agree with those expected from the PDG Fi5(1680) 
and Gi7(2190) states, assuming their mass dependencies 
are well described by the constant width Breit-Wigner 
line shape described in Eq. ( 22 ) . Allowing the = 5/2+ 



parameters to vary freely gives us better agreement and 
yields a mass around 1.95 GeV. The presence of a sec- 
ond state in the = 5/2+ wave near this mass would 
have virtually no effect on the results obtained in Fit I; 
however, it would mean that the use of a Breit-Wigner 
distribution in the energy range where both 5/2+ states 
are contributing significantly, e.g. the energy region ex- 
amined in Fit II, is invalid. 

Instead, we can employ a two-pole single channel K- 
matrix for the 5/2+ states of the form [55] 



K{s) = J2 



(24) 



Q = l 



where Wa and gapuj are the K-matrix resonance masses 
and coupling constants to the puj final state and Bi are 
the centrifugal barrier factors (see, e.g., [23]). The mass 
dependence of the amplitude is then written in terms of 
the production vector, P, and 2-body phase-space factor, 
p, as 



n^/2+{s)^p{\~ipKr 



where 



P = 



a=l 



(25) 



(26) 



with production coupling constants gap-y and 

^ ^ V(g - {wu^ + WpY){s ~ {w^ - WpY) ^27) 
s 

For this MDF, we required the = 7/2^ parame- 
ters to be within the limits quoted by the PDG for the 
Gi7(2190). One of the = 5/2+ K-matrix poles was 
required to be consistent with the Fi5(1680). The exact 
location of this pole depends on how one treats the open- 
ing of the pu) threshold (this is a single channel analysis) . 
The parameters of the second — 5 /2+ K-matrix pole 
were obtained from the MDF yielding 1930 MeV for the 
mass and 100 MeV for the width. Even with these con- 
straints, the results provide a very good description of 
the phase motion obtained from the PWE. 

In principle, poles in the T-matrix and poles in the K- 
matrix can be quite different. The relationship between 
the two can also depend on the specific K-matrix model 
employed. For these reasons, care must be taken when 
interpreting the parameters obtained for the "second" 
5/2+ resonance. While the K-matrix parameters may 
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not coincide exactly with the physical T-matrix values, 
the observed strength in this wave and its phase mo- 
tion relative to the Gi7(2190) support it having a mass 
around 1.9-2 GeV and a width of approximately 200- 
300 MeV. These values are in good agreement with the 
missing Fi5(2000) state predicted by [25]. A check us- 
ing the (unitarity violating) two Breit-Wigner prescrip- 
tion also resulted in 5/2"'' resonance parameters in this 
range. To extract precise resonance parameters for this 
state (and to confirm its existence), a coupled-channel 
analysis should be employed. 



The strengths and phase motion of the two s-channel 
waves presented in this section become noisier in the pres- 
ence of these additional waves; however, the conclusions 
drawn about resonance contributions are unaffected by 
additions to the wave set. We are unable to determine 
from our PWE's which additional waves may coincide 
with unaccounted-for physical processes; thus, we do not 
present them here. Perhaps future measurement of addi- 
tional polarization observables in this energy range might 
help determine the nature of these waves. 



2. Production helicity amplitudes 

The ratio of the helicity amplitudes for the Fi5(1680) 
was discussed in Section |IVB 2 Without employing 
a model we cannot separate out the possible missing 
Fi5(2000) production amplitudes. The ratio of the he- 
licity amplitudes for the Gi7(2190) is extracted to be 



A 



3/2 



-0.17±0.15. 



/2 



(28) 



Due to its 1-star coupling to 7p, the PDG does not quote 
a value for this ratio. 



3. Comparison to observables 

As in Fit I, we do not expect the limited number of 
waves used in Fit II to include all of the physics at these 
energies. Thus, we again do not expect to provide a per- 
fect description of the observables in this energy regime. 
Before we examine Fit II, we note that the OTL t-channel 
terms provide a good description of our forward cross 
sections; however, there are some noticeable discrepan- 
cies with the spin density matrix elements at these ener- 
gies Pl. 

Figures [9] and 10 show the differential cross sections 
and spin density matrix elements extracted from the 
PWA fits compared to our measurements TT]. As the 
energy increases, the quality of the descriptions of the 
spin density matrix elements decreases. This is partly 
due to the issues with the non-resonant model employed 
(as discussed in Section IIIB I. At these energies, the OTL 



i-channel terms begin to fail to adequately describe the 
polarization observables. The discrepancies in the back- 
wards direction could be due to the lack of inclusion of 
any u-channel terms. The effects on our results, i.e. the 
effects on conclusions drawn about resonance contribu- 
tions, due to possible issues with the non-resonant terms 
are discussed in Section IVl 

The lack of perfect description of the data signifies that 
there are other production mechanisms that are not ac- 
counted for in the fit. This is, again, expected due to the 
limited number of waves included in the PWA. Adding 
additional waves improves the description of the data. 



D. Evidence for additional resonance states 

One of the prime motivating factors in undertaking 
this study was to search for missing resonances. The 
strongest evidence for resonance contributions to uj pho- 
toproduction found in Fits I and II was for well-known 
PDG states. Suggestive evidence was also found for a 
missing i^i5(2000) state. Below we examine possible ad- 
ditional resonance contributions. 



1. Fit III: the 3/2+ wave 

Quark model calculations predict three missing reso- 
nances with = 3/2+ in the energy range of Fit I 
which couple to pu) 26J. Figure 11 shows the cross sec- 
tions and phase motion obtained if we add a = 3/2+ 



wave to the PWE in Fit I. Below 1800 MeV, the range of 
production angles over which the CLAS has acceptance 
is limited. This makes it difficult to cleanly separate 
contributions from three waves; thus, this energy range 
has been excluded from these fits. The strengths and 
phases of the — 3/2~, 5/2+ waves are virtually unal- 
tered by the addition of the extra s-channel wave. The 
cross section of the — 3/2+ wave does show some 
fairly smooth structure; however, its phase motion, rel- 
ative to the other two resonant states, is not consistent 
with a single constant width Breit-Wigner hypothesis. If 
we instead perform a MDF using the K-matrix formalism 
described in Section IV C 1 for the = 3/2+ wave, the 
phase motion between the = 3/2~, 3/2+ waves is well 
described. The poles in the K-matrix are at 1850 MeV 
and 1950 MeV. We could apply the same procedure for 
the — 5/2+, 3/2+ phase motion; however, both waves 
would have K-matrices and the number of free parame- 
ters would leave the fit under-constrained. 

Figures [12] and [13] show the comparisons of the PWA 
results with our measurements [IF. The additional 3/2+ 
wave improves the description of the spin density matrix 
elements obtained in Fit I (see Fig. J6|. It is difficult to 
make firm conclusions about the J^== 3/2+ wave. The 
observed strength suggests there is significant overlap of 
the scattering amplitude with this partial wave; however, 
this is not sufficient evidence to claim resonance contri- 
butions. The phase motion of the 3/2+ wave relative to 
the 3/2^ and 5/2+ waves is not consistent with a sin- 
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gle resonant state. This does not, however, rule out the 
existence of the multiple states predicted by the quark 
model. 



2. Limitations of the mass-independent technique 

Numerous other fits that we have performed yielded 
inconclusive evidence for states of various spin-parities. 
These fits are very similar to Fit III. Generally, smooth 
structures are found in the extracted cross sections; how- 
ever, the phase motion is inconsistent with a single res- 
onant state. It is possible that a number of resonant 
states exist which couple relatively strongly to puj. It 
is also possible that the smooth cross sections are sim- 
ply the result of overlap of various partial waves with 
unaccounted-for non-resonant terms. It would appear 
that we have reached the limits of what our technique 
can extract from our data. 

More polarization information may be required to 
cleanly extract additional resonances. Improved theo- 
retical input for the non-resonant (non s-channel) terms 
may also be necessary. While our studies have shown 
that the strong s-channel signals extracted by this anal- 
ysis are not affected by the way that the non-resonant 
terms are modeled (see Sec. [v]), this is almost certainly 
not the case for weaker signals. This is particularly true 
at higher energies, where the current theoretical models 
do a poor job of describing the new CLAS data. The 
amplitudes that are currently being generated by several 
groups from a coupled-channel approach (see, e.g., |27j ) 
may well allow for the extraction of much weaker reso- 
nance signals from these data. 

V. SYSTEMATIC STUDIES 
A. s-channel scans 

In Section |IV A[ we found that the single s- 
channel waves with the best likelihoods were 
JP = 3/2- for W < 1.85 GeV and = 5/2+ for 
1.85 GeV< W <2 GeV. We proceeded to add single 
s-channel waves to these fits to determine which had 
the best likelihoods; these were the basis for the wave 
sets used in Fits I and II. We can also examine the 
"discarded" wave sets and examine the s-channel 
contributions as a systematic check on our results. 

In the W < 2 GeV energy range, the s-channel waves 
with jP = 3/2", 5/2+, used in Fit I, had the best like- 
lihood of all two s-channel wave combinations (when 
combined with the OTL t-channel waves). In this fit, 
the contributions extracted for the two s-channel waves 
were approximately equal in size for W < 1.85 GeV. 
For all other [3/2~, J^] combinations, the extracted con- 
tribution for the — 3/2~ wave was the bigger of 
the two s-channel terms in this energy range. In the 
1.85 GeV< W <2 GeV energy range, the 5/2+ wave had 



the larger of the two s-channel contributions for every 
used for the other s-channel wave. Fits were also 
run using all two s-channel wave combinations (with the 
OTL t-channel terms) over the entire energy range. The 
contributions of the waves presented in this paper were 
consistent , regardless of which other waves they were fit 
with. 

The robustness of the results presented in this paper 
was also tested by performing the PWE's with larger 
wave sets. The presence of any additional J < 5/2 wave 
does not effect the conclusions drawn about resonance 
contributions to Fit I or Fit II. Fits run with very large 
wave sets that included all s-channel waves with J < 
9/2 also confirm the large contributions from the = 
3/2" and 5/2+ below 2 GeV and from J > 7/2 waves 
around 2.2 GeV; however, with this many waves it was 
not possible to unambiguously determine the spin-parity 
of the large J contribution. 



B. Including u-channel terms 

Another possible cause of systematic effects is our lack 
of inclusion of any u-channel terms. For W < 2 GeV, 
any w-channel contribution must be small due to the 
lack of any visible peak in the backwards cross section. 
Thus, the conclusions drawn from Fit I are independent 
of whether or not it-channel terms are included. The 
same cannot be said for the energy regime of Fit II. 
In |22j . we were able to modify the u-channel parame- 
ters of the Oh, Titov and Lee model to better describe 
our highest energy data; however, these amplitudes were 
not included in our PWA fits due to a lack of confidence 
in the assumptions used to obtain the parameters. 

We can examine what effect adding these terms would 
have on the resonance parameters extracted in Fit II. 
Figure [14] shows the phase motion obtained from Fit II 
with and without it-channel amplitudes. The agreement 
is very good in the region where both of the s-channel 
waves have strong contributions to the cross section. It 
is only in the regions where the cross section of one of the 
s-channel terms is very small that including w-channel 
terms leads to a discrepancy in the extracted phase mo- 
tion. Therefore, the conclusions drawn from Fit II re- 
garding the resonance states are unaffected by how the 
M-channel terms are modeled. 

We also note here that the likelihoods of the fits con- 
taining the u-channel terms were worse in all bins. Per- 
haps this is not surprising since the modified OTL terms 
were obtained assuming the entire backward production 
amplitude is due to u-channel mechanisms at our high- 
est energies. To obtain a better w-channel model, the 
OTL parameters should be fit including s-channel waves; 
however, to simply estimate the effects of neglecting u- 
channel terms in Fit II, these parameters are sufficient. 
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VI. CONCLUSIONS 

An event-based mass-independent partial wave analy- 
sis has been performed on data obtained using the CLAS 

at Jefferson Lab. Evidence has been found for contri- 
butions from the Fi5(1680) and 1)13(1700) nucleon res- 
onance states. These states are found to be dominant 
near threshold. The data also strongly support the pres- 
ence of the Gi7(2190) state. Suggestive evidence for an 
additional 5/2+ state with a mass around 1.9-2 GeV has 
also been found. The data shows definite strength in 
this partial wave over a very large energy range. The 
phase motion between this wave and the 017(2190) sup- 
ports the presence of a second 5/2+ state near 1.95 GeV. 
Some evidence for other states exists, although the inter- 
pretations are more difficult. The strength seen in the 
= 3/2+ wave around W = 1.8 — 2 GeV, for example, 
is not consistent with a single resonant state; however, 
we cannot rule out the possibility that multiple 3/2+ 
resonances could be contributing to our data at these 
energies. To extract additional resonance signals from 



our data, improved theoretical input for the non-resonant 

terms may be required. In particular, including the am- 
plitudes currently being generated by coupled-channel 
analyses would be highly desirable. 
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FIG. 5: (Color Online) da / dcosO'^.m.X^'b) vs cos^"„.: PWA results from Fit I (solid black line), compared to our measure- 
ments [TT]. The individual contributions from the = 3/2~ wave (dashed red line), = 5/2"*" wave (dashed-dotted blue 
line) and OTL t-channel terms (dotted black line) are also shown. The lack of data reported in the W = 1.955 GeV bin is due 
to normalization issues 1111. 
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FIG. 7: Results from Fit II: a{fib) vs W^(MeV): Total cross 
sections from all of the waves in the fit (filled squares), only t- 
channel waves (open squares), only = 5/2"*" waves (circles) 
and only = 7/2^ waves (triangles). The cross section 
extracted for = 5/2^ is consistent with the tail of a lower 
mass state (as seen in Fit I). The = 7/2~ cross section 
is indicative of a state near 2.2 GeV. The errors are purely 
statistical. 




FIG. 8: Results from Fit II: = ^67/2- ~ 4>5/2+ (radians) vs 
W(MeV): The dot-dashed line is the phase motion expected 
using constant width Breit-Wigner distributions and the pa- 
rameters quoted by the PDG for the Fi5(1680) and G'i7(2190). 
The dashed line required the — 7/2~ parameters to be 
within the PDG limits for the G'i7(2190), while allowing the 
= 5/2^ parameters to vary freely. The solid line used a 
constant width Breit-Wigner distribution for the 6*17(2190), 
but a 2-pole single channel K-matrix for the = 5/2^ wave. 
The parameters obtained from these fits are listed in the text. 
The error bars are purely statistical. 
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FIG. 9: (Color Online) da /d cos 6'^„^{iJ.b) vs cosO^.m.'- PWA results from Fit II (solid black line), compared to our measure- 
ments [TT]. The individual contributions from the = 7/2~ wave (dashed red line), = 5/2"'' wave (dashed-dotted blue 
line) and OTL t-channel terms (dotted black line) are also shown. 
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FIG. 10: (Color Online) p\[m' cosS"^.: PWA results from Fit II: poo (solid black line), Pi_i (solid red line), RepiQ (solid 
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FIG. 11: Results from Fit III: (a) (t(^6) vs ly(MeV): Total cross sections extracted from all of the waves in the fit (filled black 
squares), only t-channel waves (open squares), only 5/2"*" waves (circles), only 3/2~ waves (triangles) and only 3/2"*" waves 
(crosses), (b) A<^ = (^5/2+ ~ 03/2+ (radians) vs W(MeV). (c) = <j)3/2- ~ 03/2+ (rs^dians) vs W(MeV). The solid curves show 
the phase motion expected assuming the = 3/2"*" has Breit-Wigner parameters M3/2+ ~ 1875 MeV and r3/2+ = 150 MeV 
while locking the = 3/2~ and — 5/2'*' parameters to be those of the _Di3(1700) and -Fi5(1680), respectively. The 
phase motion obtained for the 3/2+ is not consistent with a single resonant state. The dashed line on (c) represents using 
the 1)13(1700) Breit-Wigner parameters for the 3/2~ and a single channel two-pole K-matrix for the 3/2"'". There is enough 
freedom to describe the data (see text for discussion). All error bars are purely statistical. 
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FIG. 12: (Color Online) da /d cos 6'^„^ {iJ.b) vs cos^c^. : PWA results from Fit III (solid black line), compared to our measure- 
ments [II,. The individual contributions from the = 3/2~ wave (dashed red line), = 5/2"*" wave (dashed-dotted blue 
line), the = 3/2+ wave (dashed-triple-dotted green line) and OTL t-channel terms (dotted black line) are also shown. The 
lack of data reported in the W = 1.955 GeV bin is due to normalization issues |1I] . 
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FIG. 13: (Color Online) p\, j^ji vs cosS" PWA results from Fit III: poo (solid black line), pi_i (solid red line), RepiQ (solid 
blue line), compared to our measurements [TT]. No fit was performed in the W — 1.955 GeV bin, see text for details. 
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bars on both phase motion plots are purely statistical. 



